Nonlinear Breathing-like Localized Modes in Cgo Nanocrystals 
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We study the dynamics of nanocrystals composed of Ceo fuUerene molecules. We demonstrate 
that such structures can support long-lived strongly localized nonlinear oscillatory modes, which 
resemble discrete breathers in simple lattices. We reveal that at room temperatures the lifetime of 
such nonlinear localized modes may exceed tens of picoseconds; this suggests that Ceo nanoclusters 
should demonstrate anomalously slow thermal relaxation when the temperature gradient decays in 
accord to a power law, thus violating the Cattaneo-Vernotte law of thermal conductivity. 



I. INTRODUCTION 

The fuUerene era was started in 1985 with the dis- 
covery of a stable Cgo cluster and its interpretation as 
a cage structure with a soccer ball shape [1]. By now, 
many structural, electronic, and vibrational properties 
of fuUerenes have been studied in detail. In spite of 
the rapidly increasing interest in new forms of fullerenes, 
icosahedral Cgo, the "most beautiful molecule" [2], re- 
mains at the focus of active research as a prototype 
fuUerene system. FuUerene molecule Ceo keeps a promi- 
nent position primarily for the high symmetry which 
greatly facilitates theoretical computations and the in- 
terpretation of experimental data, but also because it is 
only fuUerene which can be grown in the form of large and 
nearly perfect single crystals. At room temperature, the 
diffraction experiments [3] show that the Cgo molecules 
are centered on sites of a face-centered-cubic (fee) Bra- 
vais lattice. Each Cgo molecule in the solid-state phase 
is coupled to other Cgo molecules through weak Van-der 
Waals forces. 

By now, the internal oscillations of the fuUerene 
molecule Cgo are well studied, for the case of an isolated 
molecule and also for Cgo crystalline structures. Various 
theoretical and numerical methods have been developed 
and employed. However, the dynamics of fullerenes in the 
nonlinear regime and a role of nonlinear localized modes 
in such structures remain largely unexplored. 

The main aim of this paper is to study strongly lo- 
calized nonlinear oscillatory modes in such nanocrystals 
which resemble the discrete breathers of much simpler 
nonlinear lattices. 

High symmetry of Cgo molecules allows the existence 
of internal nonlinear oscillations which do not alter sub- 
stantially the spherical shape of the fuUerene molecule. 
Such oscillations are excited by rotational degrees of 
freedom of pentagons and hexagons of the structure, 
with a variation of valent bonds. Resonant coupling 
of such high-amplitude oscillations with the neighboring 
Cgo molecules will remain weak, because due to conser- 
vation of a spherical shape the molecules will interact 
as large particle through Lennard- Jones potentials. This 
observation suggests that such structures can provide an 
ideal environment for the existence of nonlinear modes 



strongly localized in the space [4, 5]. 

Here we demonstrate that the Cgo crystals can sup- 
port long-lived strongly localized nonlinear modes which 
closely resemble the so-called discrete breathers in sim- 
ple nonlinear lattices, and affect dramatically the heat 
relaxation of composite fuUerene structures that obeys a 
power low instead of an exponential law. 



II. MODEL AND SIMULATION METHODS 

To analyze the dynamics of of localized nonlinear 
modes in three-dimensional Cgo crystals, we consider a 
cubic nanocluster of Cgo molecules, as shown in Fig. 1 (a). 

Inside of each Cgo molecule of the nanocluster two 
neighboring carbon atoms create a valent bond. We de- 
scribe the energy associated with deformation of the va- 
lent bond created by two carbon atoms by the following 
interaction potential 



Ui{p) = ei{exp[-ao(/5 - Po)] - 1}^ 



(1) 



where p is the length of the valent bond, ei — 4.9632 eV 
and Po — 1.418 A is the energy and the equilibrium 
length of the bond, parameter ao = 1.7889 A~^. Each 
carbon atom is placed at the vertex of three planar valent 
angles. The corresponding deformation energy of a plane 
valent angle created by three atoms can be described by 
the interaction potential of the form. 



t/2((/5) = 62 (cos V5+ 1/2)2 



(2) 



where ip is value of the valent angle (equilibrium value 
of the angle ipo — 2n/3), energy £2 = 1-3143 eV. Each 
valent bond is simultaneously belong to two planes, so 
that the deformation energy of the angle formed by two 
such planes can be described by the following interaction 
potential. 



(1 



(3) 



where 4> ~is the corresponding angle (in equilibrium, (jj = 
0), and the index z = 3, 4, 5 describe the type of the 
valent angle. Energy £3 = £4 = 0.499 eV, £5 <^ £4, so 
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FIG. 1: (a) General view of a cubic nanocluster composed of 5 x 5 x 5 fullerene molecules Ceo (buckyballs C60-Ih), (b) isolated 
molecule Ceo and (c) dimensionless density of its frequency spectrum for thermal vibrations at T = 3K. 



that the latter contribution to the total energy can be 
neglected. 

More detailed discussion and motivation of our choice 
of the interaction potentials Ui,...,Ui can be found in 
Ref. [6]. Such potentials have been employed for mod- 
eling of thermal conductivity of carbon nanotubes [7, 8] 
graphene nanoribbons [6, 10] and also in the analysis of 
their oscillatory modes [9, 11, 12]. 

We describe the interaction of carbon atoms belong- 
ing to different Ceo molecules by the Lennard-Jones 12-6 
potential 

ULjir)^4eU{'r/rr-{a/r)% (4) 

where r - distance between the atom centers, energy 
eLj = 0.004556 cV, and parameter a = 3.851 A. 

To verify our results, we also employ the simpler model 
where all interacting potentials between the atoms are 
described, instead of the set of the potentials Ui,..., and 
C/4, by the Brenner potential [13]. 

III. INTRAMOLECULAR MODES 

The high symmetry of an isolated, isotopically 
pure Cgo molecule in its equilibrium configuration 
imposes very strong constraints on the form of the 
normal mode displacement patterns. All 60 sites 
are symmetry-equivalent, and the point group is the 
full icosahcdral group Iji. Although there are 180 



degrees of freedom (3x60) for each Cgo molecule, 
the icosahedral symmetry gives rise to a number of 
degenerate modes, so that only 46 distinct mode 
frequencies are expected for the Ceo molecule. All 
vibrational modes of Ceo arc distributed [14, 15] as 
2Ag+A„+4Tig+5Ti„+4T3g+5T3„+6Gg+6G„+8Hg+7H, 

Dynamics of an isolated Ceo fullerene molecule is de- 
scribed by a system of 3 x 60 nonlinear equations of the 
standard form, Afii ~ —dH/du, where M = 12 ■ 1.6603 ■ 
10^^^ kg - mass of the carbon atom, system Hamilto- 
nian is = A/(u, u)/2 + E{u), and 180-dimensional 
vector u — {{xn,yn, Zn)}^=i defines the coordinates of all 
atoms of the complex molecule. Potential energy of the 
molecule E{u.) is composed by all types of interactions 
described by the potentials Ui,...,U4 discussed above. 

To find the equilibrium stationary state of the Ceo 
molecule, we need to solve numerically the minimization 
problem for the energy functional, E{u) — > min. In the 
model under consideration, the molecule Ceo is described 
by the set of coordinates Uq , and it has a shape of perfect 
icosahedral [see Fig. 1 (b)] with the radius R = 3.514 A, 
for which all valent bonds have the equilibrium length 
P = Po- 

Linearizing the model near the equilibrium station- 
ary state, we obtain a system f linear equations Alii = 
Au, where A is 180 x 180 dimensional matrix, A = 
{dE /duiduj\u=uo}i^i^fJ^i- Diagonalizing the matrix A 
we find all 180 eigenmodes of the molecule. The first 
six zero- frequency modes correspond to rigid translations 
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TABLE I: Values ijj (in cm ^) for the 46 distinct normal mode frequencies of an isolated Ceo molecule (n is mode degeneracy). 
Notations of the modes are from Ref. [15] . 
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and rotations of the molecule. Other 174 modes are gen- 
erated and, in view of their degeneracy, there are exists 
46 distinct modes - see Table L 

The frequencies of all oscillatory eigenmodes of an 
isolated Cgo molecule arc summarized in Table I. All 
frequency are located in the interval 162.1 < io < 
1579.3 cm~^. Definitely, the model with only 6 param- 
eters is not able to predict exactly the frequencies ob- 
tained in the framework of the first principles calcula- 
tions or experimental results, see [15]. Nevertheless, this 
model demonstrates a correct structure of the eigenmode 
frequencies and predict rather well the frequency band 
measured experimentally, 272 < w < 1575 cm^^. 

It is important to mention that the use of the Bren- 
ner potentials [13] instead of the set of the potentials 
Ui,...,U4 leads much broader frequency band, 182.9 < 
u! < 1738.9 cm~^. In addition, the Brenner poten- 
tial do not predict correctly the degeneracy of the fre- 
quency spectrum, so that the equilibrium shape of the 
Ceo molecule deviates from the icosahedra shape. 

In order to model thermal oscillations of the molecule, 
we start from its equilibrium state and excite each of its 
eigenmode with the energy ksT. The frequency spec- 
trum of the Ceo molecule at the temperature T = 3K is 
shown in the Fig. 1 (c). As expected, the spectrum is dis- 
crete, and the number of peaks coincides with the num- 
ber of eigenmodes, whereas the peak positions coincide 
with the corresponding eigenfrequencies. At low temper- 
atures, the oscillations are linear, so that the spectrum 
corresponds to a system of uncoupled linear oscillators. 

If we increase the temperature to T = 30K, the discrete 
nature of the spectrum remains unchanged, but high- 
frequency modes demonstrate some energy exchange and 
interaction. For much higher temperatures such as T = 
300K the discrete nature of the spectrum is lost and the 
oscillations become strongly nonlinear demonstrating in- 
teraction with strong energy exchange. 



IV. LOCALIZED NONLINEAR MODES 

The Ceo molecules occupy fee sites in the solid phase, 
and at room temperature, the rapid rotational diffusion 
of each molecule leads to an effective fee crystal structure. 
To study localized modes of the fullerene molecule placed 
in a crystal, we consider a cubic nanocluster composed 
of 5 X 5 X 5 molecules Ceo, as shown in Fig. 1 (a). 

First, we find the ground state of the nanocluster by 
minimizing its energy. Next, for the central molecule we 
excite its n-th eigenmode {n ~ 7, ...,180) with the en- 
ergy Eq and study the decay of the kinetic energy in 
the cluster. Due to the nonvalent coupling between the 
fullerenes, the energy is transferred from one molecule to 
another. To model the energy spreading in an effectively 
infinite crystal, we introduce lossy boundary conditions, 
when all surface molecules experience damping with the 
relaxation time tr = 20 ps. If the excitation of the cer- 
tain eigenmode lead to the creation of a localized state, 
or breathing-like mode, the kinetic energy will not vanish 
but instead will approach a certain nonzero value. Oth- 
erwise, the initial excitation will vanish completely, so 
that the kinetic energy of the central molecule Ek — !• 
when t — > oo. 

First, we excite with the energy Eq = 0.2 the n- 
th mode {n = 7, ...,180) of the central molecule of the 
fullerene crystal. Depending on the type of the modes, 
our numerical results reveal three different scenarios of 
the mode relaxation, as shown in Fig. 2 (a-c). The first 
scenario is the rapid energy relaxation and decay of the 
mode oscillations (see curves 1, 4, and 7); the second sce- 
nario is slow relaxation (curves 2, 5, and 8), and the third 
scenario is the generation of stationary spatially localized 
nonlinear modes (curves 3, 6, and 9). 

The first and second scenarios are observed for the 
modes associated with a periodic change of the spherical 
shape or the size of the Ceo molecule. The most interest- 
ing third scenario is observed for the modes associated 
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FIG. 2: Evolution of the kinetic energy Ek of the central 
fullerene Ceo molecule in nanocluster composed of 5 x 5 x 5 
molecules packed in the form of a fee cubic crystal. At the 
initial moment t = only one eigenmode is excited for the 
central molecule with the energy Eo = 0.2eV. (a) curve 1 show 
the decay of the mode Hg(l) (frequency uj — 162.2 cm~^), 
curve 2 - the mode Hu(4) [uj = 501.8 cm~^), and the curve 
3 - the mode T3g{l){iJ = 442.3 cm"^); (b) curve 4 gives the 
dependence for the mode Ag{l){uj = 418.6 cm~^), curve 5 - 
the mode T35(4) (a; = 1267.6 cm~^), curve 6 ~ the mode 
Ti9(2) (w = 675.7 cm-i); (c) curve 7 describes the decay of 
the mode A^ (frequency u — 776.0 cm~^), curve 8 - the mode 
09(5) {uj = 1437.5 cm"^), and the curve 9 - the mode Tig(3) 
{lu = 1454.9 cm"^). 



with the motion of atoms along the spherical surface 
so that the molecule keeps its form almost unchanged. 
These are the eigenmodes T3g(l), Tig(2), and Tig(3) 
presented in Fig. 3. For all those modes, the molecule 
conserves its spherical shape and therefore the resonant 
interatomic interaction is reduced. 

Figure 4 shows the spatial distribution of the kinetic 
energy of the nonlinear localized state generated by the 
excitation of the eigenmode T3g(l). The main part of 
the energy is concentrated at the central Ceo molecule, so 
this state resembles discrete breathers studied earlier for 
simpler lattices. The nonlinear states generated by the 
eigenmodes Tig (2) Tig (3) are localized even stronger. 
For E > 0.01 eV, all such modes represent strongly an- 
harmonic oscillations with the frequency decaying with 
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FIG. 3: Selected eigenmodes of an isolated Ceo molecule 
which lead to the generation of spatially localized nonlinear 
modes (breathers) of the fullerene crystal, (a) mode T3g(l) 
(frequency u = 442.3 cm^^); (b) mode Tig (2) (frequency oj = 
675.7 cm"^); (c) mode Tig(3) (frequency lu = 1454.9 cm"^). 



(a) . " . (b) 




FIG. 4: Spatial profile of the kinetic energy E of 
the breathing-like nonlinear localized mode of the central 
fullerene molecule in the nanocluster of 5 x 5 x 5 of Ceo 
molecules in the planes (a) xy and (b) xz, intersecting the 
central molecule (shown by a black circle) For the generation 
of such breathers, the eigenmode T3g(l) was excited with the 
energy E = 0.2 eV. 

the mode amplitude. The mode Tig (2) carries the largest 
amount of energy that can reach 1 eV. 

It is important to mention that the model employing 
the Brenner potentials predicts the existence of a larger 
number of localized modes. 

In a thermalized nanocrystal all localized modes have a 
finite life time. Their decay occurs mainly due to the in- 
teraction with other modes of the same molecule. There- 
fore, in order to define the lifetime of a nonlinear mode, it 
would be sufficient to study the dynamics of an isolated 
Cgo molecule with one strongly excited mode correspond- 
ing to the nonlinear localized state and all other modes 
weakly thermalized. 

In Fig. 5 we show the interaction of one large- 
amplitude eigenmode T3g(l) with all other weakly ther- 
malized eigenmodes of the fullerene Cgo molecule. As fol- 
lows from those results, the energy of the localized mode 
decays monotonically exciting other modes. When the 
initial energy of the T3g(l) mode is £^1 =0.2 eV and the 
energy of all other modes is Eq — 0.002 eV, the localized 
mode loses a half of its energy for to = 1.4 ns, for Eo = 
0.004 eV "this time is to = 0.8 ns, for Eo = 0.008 eV - 
to = 0.4 ns, and for Eo = 0.016 eY ~ to = 0.1 ns. Highly- 
excited mode Tig (2) with the energy Ei = 0.2 eV losses a 
half of its energy for to = 0.8 ns if energy Eo = 0.002 eV, 
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FIG. 5: Interaction of a single large-amplitude mode 
T3g(l) with all other weakly thermalized eigenmodes of Ceo 
molecules. As the initial condition ai t = we excite the 
eigenmode T3g(l) with the frequency u = 442.3 cm~^ and 
the energy Ei — 0.2eV, and all other eigenmodes with the 
energy Eq = 0.002eV. Shown is the evolution of the energy 
distribution over the frequencies. The peak corresponds to a 
large-amplitude mode (shown in red). 

for to = 0.5 ns if Eq = 0.004 eV, for to = 0.3 ns if 
Eo = 0.008 eV, and for to = 0.1 ns if Eq = 0.016 cV. 
The smallest lifetime is shown by the eigenmode Tig (3), 
which for the input energy Ei = 0.2 eV and Eo = 0.002 
looses a half of its energy for to ~ 0.1 ns. 

Therefore, our numerical results demonstrate that in 
a crystal formed by Ceo molecules the spatially localized 
nonlinear modes are long-lived excitations, and at room 
temperature their life time is of the order of several pi- 
coseconds. 



V. THERMAL RELAXATION OF Cgo CRYSTAL 

The existence of long-lived nonlinear localized collec- 
tive modes in the nanocluster of Cgo fullerene molecules 
should affect substantially their physical properties. If 
the thermal energy can be concentrated for long time in 
some localized spots, such a crystal is expected to demon- 
strate low thermal conductivity. To confirm this idea, we 
consider thermal relaxation of a cubic nanocluster com- 
posed of 10 X 12 X 12 Ceo molecules, as shown in Fig. 6. 

The 180-dimensional vector Unmk define the coordi- 
nates of the atoms of (n, m, fc)-th molecule of the nan- 
ocluster {n = 1,...,10; TO = 1,...,12; k = 1,...,12). 
For thcrmalizing the nanocluster, we solve the system 
of Langcvin equations 

Milnmk = -dH/dUn,nk - rA/u„,„fe + E^mk, (5) 

where F = l/t,. is the Langcvin collision frequency with 
tr ~ 0.1 ps being the corresponding particle relaxation 
time, and S„,„fc = {Cnmfe,i}ili is a 180-dimensional 
vector of Gaussian distributed stochastic forces de- 
scribing the interaction of atom with the thermostat 
with correlation functions (^„imifci,i(ii)^„2m2A;2,j(t2)) = 




FIG. 6: Relaxation the thermal gradient AT between the op- 
posite sides of the cubic cluster of 10 x 12 x 12 fullerene Ceo 
molecules with time t Left insert shows this dependence in the 
log-scale, with a straight line corresponding to a power-law de- 
pendence, AT '--^ t~^. Right insert shows a cubic nanocluster 
with an arrow indicating the heat flow. Dependencies 1 and 
4 correspond to the approximation of rigid molecules when 
their internal dynamics is neglected, whereas curves 2 and 5 
show the results of the full molecular dynamics simulations. 
Curves 3 and 6 show the results for the model employing the 
Brenner potentials. 



Sij6{t2 — ti), Tnmk is the 
temperature of the molecule with the site (n, to, k). 

To obtain the initial thermalized state of the nanoclus- 
ter we solve numerically the Langevin equations (5) with 
the initial conditions corresponding to the equilibrium 
state. We run the simulations for t = lOt,., and then 
switch off the coupling with the thermostat, i.e. neglect- 
ing the last two terms in Eqs. (5), in order to study the 
thermal relaxation of the nanocluster. 

We thermalize the nanocluster asymmetrically, to have 
only one side heated to the temperature T = 500 K, 
whereas keeping all other atoms at the temperature T = 
300K. To achieve this, we solve the system of Langevin 
equations at Tnmk = 500 K for the atoms with the index 
k = 1 but for all other atoms with the indices k = 2, ...,12 
we put Tnmk = 300 K. In this case, we achieve a tem- 
perature different of AT = 200 K between the opposite 
sides of the Ceo nanocube. As the next step, we study 
the thermal relaxation of the temperature gradient. 

As a result of thermal energy exchange between the 
fullerene molecules, temperature difference should decay 
monotonically. The crystal of fullerene molecules Ceo 
has a three-dimensional structure so it should obey the 
Fourier law for thermal conductivity. Because time docs 
not enter the Fourier law, it does not allow to deter- 
mine the temporal relaxation of the thermal gradient, 
for which we should employ the well-know generaliza- 
tion known as Cattaneo-Vernotte (CV) law [16, 17]. In 
its one-dimensional version, this law can be formulated 
through the differential equations, 

(l + r|),^^-.VT, (6) 
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where k is the standard coefficient of the thermal con- 
ductivity and T is the characteristic relaxation time of 
the system. In according to Eq. (6), the temperature 
gradient AT should decay in accord with an exponential 
dependence. 

Dependence of the temperature gradient AT vs. t is 
shown in Fig. 6. As follows from this dependence, the 
temperature difference decays very slow, in accord with 
the power-law dependence, so that it does not obey the 
Cattaneo-Vernotte law. In order to understand the rea- 
son for such an anomalous thermal relaxation as well as 
the role played by the internal dynamics of the molecules, 
we replace the Ceo molecules by rigid balls and solve the 
same problem assuming that the molecules in the crys- 
tal interact between each other as structureless particles. 
Our numerical results show that in this approximation 
the relaxation of the temperature difference AT follows 
an exponential dependence, in accord with the Cattaneo- 
Vernotte law. Therefore, the anomalously slow relaxation 
can be attributed to the excitation of long-lived nonlinear 
oscillations of the Cgo nanocluster. 

VI. CONCLUSIONS 

We have applied the molecular-dynamics numerical 
simulations to study the nonlinear dynamics and ther- 



mal relaxation of crystalline structures composed of Cgo 
fuUerene molecules. We have revealed that such complex 
nanoclusters support a special type of strongly localized 
nonlinear modes resembling discrete breathers of simple 
nonlinear lattices. For such modes, the kinetic energy 
is localized primarily in the rotational modes of a sin- 
gle Cgo molecule and it decays slowly inside the fullerene 
nanocrystal. Importantly, our results have demonstrated 
that at room temperatures the lifetime of such nonlin- 
ear localized modes may exceed tens of picoseconds. The 
existence of such long-lived localized nonlinear states in 
the nanoclusters explains the anomalously slow thermal 
relaxation observed in such structures when the tempera- 
ture gradient decays in accord with the power but not ex- 
ponential law, thus violating the Cattaneo-Vernotte law 
of thermal conductivity. 
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